SOME CONSIDERATIONS ON THE DERIVATION 
OF THE NONLINEAR QUANTUM BOLTZMANN EQUATION. 



D. Benedetto*, F. Castella+, R. Esposito # and M. Pulvirenti* 

Abstract. In this paper we analyze a system of N identical quantum particles in a weak- 
coupling regime. The time evolution of the Wigner transform of the one-particle reduced 
density matrix is represented by means of a perturbative series. The expansion is obtained 
upon iterating the Duhamel formula. For short times, we rigorously prove that a subseries of 
the latter, converges to the solution of the Boltzmann equation which is physically relevant 
in the context. In particular, we recover the transition rate as it is predicted by Fermi's 
Golden Rule. However, we are not able to prove that the quantity neglected while retaining 
a subseries of the complete original perturbative expansion, indeed vanishes in the limit: 
we only give plausibility arguments in this direction. The present study holds in any space 
dimension d > 2. 



1. Introduction. 

A large quantum particle system in a rarefaction regime should be described by a 
Boltzmann equation. However, while the rigorous validity of the Boltzmann equation has 
been proved for classical systems for short times [L] or globally in time for special situations 
[IP] (see Ref. [CIP] for further comments), there is no rigorous analysis for the equivalent 
quantum systems. 

The problem is physically relevant because quantum effects, although usually negligible 
at ordinary temperatures (except for few light molecules), happen to play a role in the 
applications at mesoscopic level. We refer, for example, to the treatment of electron 
gases in semiconductors (for physical references, see the textbooks [RV], [AM], [Ch], as 
well as [Bo], [CTDL] - see also the articles [Fi] or [Co] - see [MRS] for a mathematically 
oriented presentation). Therefore, establishing a well founded quantum kinetic theory is 
certainly interesting not only from a conceptual viewpoint but also from a practical one. 
In fact, kinetic descriptions for quantum systems, beside dilute gases, include dense weakly 
interacting systems, as e.g. the electron gas in semiconductors, whose classical analogues 
rather yield diffusion processes. 

One pragmatic way to introduce the Quantum Boltzmann equation (see e.g. [CC]) is to 
solve the scattering problem in Quantum Mechanics and then to replace, in the classical 
Boltzmann equation, the classical cross section with the Quantum one. 

A better logically founded approach is to derive an evolution equation for the Wigner 
transform of a quantum state associated to a dilute particle system. Working on this 
equation, one can hope to recover, at the quantum level, the same physical arguments 
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than those used at the classical level to obtain propagation of chaos and a suitable kinetic 
description for the one particle distribution function. This is the strategy we adopt in the 
present paper to treat quantum iV-particle systems. We refer to the textbook [CIP], or 
the article [L] for the analysis of the classical case. 

Summarizing, in the present paper, we consider a quantum iV-particle system and 
represent the time evolution of the one-particle Wigner function, when N — > oo in the 
weak-coupling limit, in terms of a perturbative series expansion. On the basis of some 
heuristic arguments developed in the next section, we neglect some terms and consider only 
a subseries which is proven to converge, for short times, to the solution of the Boltzmann 
equation with a suitable cross section computed by quantum rules (namely the Fermi 
Golden Rule) . Therefore, the present analysis is not a rigorous derivation of the Quantum 
Boltzmann equation, but we hope it constitutes a step in this direction. Our main result 
is Theorem 3.1 below. Although we work in dimension 3 the present statements are easily 
extended in any space dimension d > 2. Our analysis heavily relies on stationary phase 
computations (see Proposition 3.3), as well as appropriate representations of the various 
solutions of the hierarchies we need to handle. We leave further comments to the last 
section and conclude the present one by establishing the model and the scaling. 

We consider a iV-particle quantum system in M 3 . We assume the mass of the particles, 
as well as h, to be one. The interaction is described by a two-body potential 4> so that the 
potential energy is: 

U(x 1 . ..x N ) = y^0(ffi - Xj). (1.1) 

i<j 

The Schrodinger equation reads: 

id t V(X N ,t) = ~A N ^(X N ,t) + U(X N )^(X N ,t) (1.2) 

where A^r = X^ili Aj is the Laplacian with respect to the Xi variables, and X^ is a 
shorthand notation for x\ . . . xn- 

We rescale the equation according to the hyperbolic space-time scaling 

x — > ex , t — > et (1-3) 

and simultaneously we rescale also the potential <j) — > y/e<j). Hence the resulting equation 
is 

e 2 

ied t y £ (X N , t) = --A N V e (X N , t) + U £ (X N )-$ £ (X N , t), (1.4) 

where 

U £ (x 1 . . . X N ) = ^ 4> e {Xj - Xj) (1.5) 

and 

^ = (1.6) 

e 

Note that \I/ e (A^v,t) is fully determined by Eq. (1.4) and the initial datum which will be 
specified later on. 
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We want to analyze the limit e — > in the above equations, while keeping 

N = e~ 3 . (1.7) 

This kind of limit is usually called weak-coupling limit. Another possible scaling to be 
considered is the low-density limit. In this case cj> is unsealed but N = 0(e~ 2 ). In the 
classical context this is nothing but the Boltzmann-Grad limit (see e.g. [CIP]). In the 
present paper we will only be concerned with the weak coupling limit which is, to some 
extent, technically easier. 

We now introduce the Wigner function: 

W N (X N , V N ) = (J^j j dY N e ^-^r (X N + £ -Y N W{X N - £ -Y N ). (1.8) 
A standard computation yields: 

(d t + V N - W N )W N (X N ,V N ) = -^{T £ N W N )(X N ,V N ) (1.9) 

where Vn • Vat = J2iLi v i ' V Xi and (d t + Vn • Vat) is the usual free stream operator. Also, 



we have introduced 



with 



(T N W N )(X N ,V N ) = Yl ( T k,eW N )(X N ,V N ), (1.10) 

0<fc<^<Af 



m A W N )(X N ,V N ) = ](^) 3N J dY N J dV N j Y "< v "- y ^W N {X N ,V N ) 

x k -x e 1 , A f x k - Xf 1 

- 7;KVk ~ Vi) - <f>\ + o iVk - Vi) 



e 2 W ~ ) T \ e 2 
In other words, 

(n A w N )(x N ,v N ) = -i y: ° [ ^kh)^-^ 

*=±* J w (i.ii) 

„.«, oh ah . 

W (x 1 ,V 1 ,...,X k ,V k + — ,...,x N ,v N ). 

The operator describes the "collision" of particle with particle £, and the total 

operator T N takes all possible "collisions" into account. Here and below, / denotes the 
Fourier transform of /, normalized as follows: 

f(h) = {T x f){h) = [ dxf(x)e- ih - x , (1.12) 
/(*) = I (1.13) 
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We now introduce the Wigner transform of the partial traces according to the formula, 
for j = 1. ...,N- 1: 

fj f (Xj,V j ) = J dx j+1 ...dx N J dv j+1 ...dv N W N (Xj, x j+1 . . . x N ; Vj, v j+1 . . . v N ) 

(1.14) 

Obviously, we set f$ = W N . 

From now on we shall suppose that, due to the fact that the particles are identical, the 
objects which we have introduced (\I/ e , W N , fj*) are all symmetric in the exchange of 
particles. 

Proceeding as in the derivation of the BBKGY hierarchy for classical systems (see [CTP]), 
we readily arrive at the following hierarchy of equations (for 1 < j < N): 

(d t + pv k . v fe )/f = -LTjfj + *L-Lq +lf » +1 . (Lis) 

The operator C| +1 is defined as: 



and 



C e kJ+1 f j+1 (x 1 ...x J ;v 1 ...v j ) = -i a J J dx j+1 J dv j+1 <j>(h) 



(1.17) 



The operator C| - +1 describes the "collision" of particle k, belonging to the j-particle 
subsystem, with a generic particle outside the subsystem, conventionally denoted by the 
number j + 1 (this numbering uses the fact that all the particles are identical). The total 
operator C| +1 takes into account all such collisions. As usual ([CTP]), equation (1.15) 
shows that the dynamics of the j-particle subsystem is governed by three effects: the free- 
stream operator, the collisions "inside" the subsystem (the T term), and the collisions with 
particles "outside" the subsystem (the C term). 

We fix the initial value {fj}f =1 of the solution {/^(t)}^ and we assume for simplicity 
that {ffjjLx is factorized, that is, for all j = 1, N 

f? = /o® \ (1-18) 



where fo is a one-particle Wigner function which we assume also to be a probability 
distribution. We remind that the quantum state, whose Wigner transform is a general 
positive fo, is not in general a wave function but rather a density matrix. As a consequence 
the evolution equation we have to use is not the Schrodinger equation (1.2) but rather 
the Heisenberg equation for the density matrix. In both cases the corresponding Wigner 
equation is (1.9). 
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One can try to handle the hierarchy (1.15) as for the Boltzmann-Grad limit for classical 
systems, namely to study the asymptotic behavior of the solution expressed in terms of 
the series expansion for 1 < j < N, obtained upon iterating the Duhamel formula, 

Sint(tl ~ tz)Cj + 2 ■ ■ ■ Sf nt (t n -i — tn)Cj +n Sf nt (t n ) f 'j +n . 

Here Sf nt (t)fj is the j-particle interacting flow, namely the solution to the initial value 
problem: 

(d t + V 3 ■ V 3 )Sf nt (t)f 3 = ^TjSlrtWfj, 

If we expand Sf nt {t) as a perturbation of the free flow S(t) defined as 

(S(t)f J )(X J ,V J )=f J (X J -V J t,V J ), (1.21) 

we find 



1 / f Tl PTm-l 



(1.22) 



S(n - r 2 m . . . 5(r m _! - r m )T e 5(r m )/ J , 



Inserting (1.22) into (1.19), the resulting series contains a huge number of terms. However, 
we claim that many of these contributions are negligible in the limit. In the next section, 
we shall give heuristic arguments in that direction. We first make our point more precise. 
We write shortly (1.19) and (1.22) as 

N-j 

fN \ r qe /~ie qe s~ie qe (~ie qe 

Jj — 2.^1 °int K - y j + l°int K - y j+2 ' ' ' ° int^ j +n°int > 

n=0 



m>0 

with an obvious abuse of notation. Here the tilde above C and T terms absorbs the 
normalization by N/ y/e in (1.19), respectively in (1.22). Now, the insertion of (1.19) 

into (1.22) readily gives 

/;=EE E • • • E (sf^sq^sf^r^sq^ ... 

n=0m >0mi>0 m„>0 {l.ZS) 

■ ■ ■ (Sfj+^r^SC^iSfj+^S . 

In view of the expansion (1.23) of f^, we first claim that all the relevant terms in (1.23) 
is that corresponding to 

mo = , mi = m2 = ■ ■ ■ = m n = 1 . 
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Hence as e — > 0, is asymptotic to, 

N-j 

~ E SCj +1 ST j+1 SCj +2 STj + 2 ■ ■ ■ STj +n _ 1 SCj +n ST j+n S . (1-24) 

n=0 

On the more, expanding each "collision" term in (1.24) into Cj +1 =^2 3 r=1 C^j +1 , and TJ = 
Yjl i=i^r ii we a ^ so c l a i m that is in fact asymptotic to, 

N-j j j+1 j+n—1 

if ~ ^ ^ ^rij + l^rij + l X] 5<C 'r 2 ,j+2 ^2,^+2 • • • E SC* n j +n ST rn j +n S . 
n=0 n = l r2 = l r n = l 

(1.25) 

In other terms, we claim that the dynamics of the j-particle subsystem is only made up of 
collision/recollision events in the asymptotics e — > 0. 

Summarizing the above claims, we define the sequence {f^jLi by 

~[ W £ ) Jo Jo Jo Jo 

t„-i /•*„ J J + 1 J+n-l 

dt„ / rfr n E E • • • E - ti)C £ ruJ+1 S( tl - n)T^ j+1 (1-26) 

" / ° r 1 = l r 2 = l r„ = l 

5(n - t 2 )C e r2 j+2 S(t2 - r 2 )T r e i J - +1 . . . 

' ' ' S(r n -i — t n )Cr n j +n S(t n — T n )T^ n j +n S(T n )fj +n . 

Obviously is a subseries of the series expansion (1.23) defining the true value f^. It 
relates the value of the right-hand-side of (1.25). The remainder part of this paper is 
dedicated to the rigorous proof that fj(t) converge to (f(t))® J for short times, where /(£) 
solves the Boltzmann equation with the suitable cross-section. The result will be precisely 
established and proved in Section 3. 

Note that unfortunately, we are not able to rigorously prove the asymptotics (1.25): neither 
can we prove reasonable uniform bounds on the relevant series expansions, nor can we even 
prove that the asymptotics (1.23) holds term- by-term . 

2. Heuristic considerations. 

In this section, we give some reasons justifying the claimed asymptotics (1.24) and 
(1.25). These are illustrated upon analyzing the lower order terms in the true expansion 
(1.23). This section thus heuristically justifies the fact that we restrict ourselves with the 
mere analysis of in the present paper. 

Also, we analyze the term corresponding to n = 1 in the expansion (1.25) defining ff. 
We prepare in this way the analysis of the complete series expansion performed in the next 
section. 
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To have an idea of the various orders of magnitude of the terms appearing in (1.19) 
once Sf nt is expanded (see (1.23)), we consider the first two terms in (1.19) expanding Sf nt 
up to the first order. As a result we have the following five terms which we are going to 
analyze: 

X = S(t)ff , (2.1) 
X x = ^d- J dh S(t - f i)CJ +1 S(f , (2-2) 

X 2 = 4= f dn S(t - r 1 )T/5(r 1 )/ J , (2.3) 
Jo 

X 3 = —j 1 J dn J dh S(t - r 1 )T/5(r 1 - t^C^S^f^ , (2.4) 

l 4 = —y^ j dh j d n S(t - h)q +1 S{t! - n)T' i+1 S( Tl )f° +1 



E E V- 

r=l l<s<e<j+l 



(2.5) 



t rti 



I rAs = NJi dh j d^sit-tJC^Sih-rjTlsSin)/^. (2.6) 



where 

q-r,i,s ^ 

As we shall see, the terms Xj, i = 1, 2, 3 are negligible in the limit e — > 0. This illustrates 
(on a particular example) the claim (1.24) above. Also, all the contributions to X4 but 
that for r = I and s = j + 1, corresponding to a collision/recollision event, are equally 
vanishing. This illustrates the claim (1.25) above. 

Clearly Xq does not require any asymptotic analysis. 
For Xi we have: 

= ^rE E ° [d tl f d Xj+ , j dv i+ , f ^ (fe )e<M--- <•->-.,« ) 

n , cx/i cr/i 
/ J+ i(xi - wit, . . . , x r - iv* + -y^i) • • • > - ^i+i^i 2~*i) 

ah 

V 1 ,...,V r - Vj+i) . 

(2-7) 

Now, upon using the stationary phase Theorem for the variables a^+i and h in formula 
(2.7), we readily obtain the following asymptotics, valid if f® +1 and 4> are smooth and 
decaying enough (which we assume), 

TV" if ^ f^f 

X x = — ~^J- x zV x V V a / efti / dv i+ i0(O) 

V £ V r=la=±l J ° J 

/jVl^l _ • • • - frt, • • • ,X r - 1> r (t - tl) - fj + ltl, (2.8) 

v 1 ,...,v r ,...,v j+1 ) +0(e 4 ) . 
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In other words, the integral in (2.7) concentrates on the set (£-£i), h = 0, 

as e — > 0, and it has size £ 3 + 0(e 4 ). Recall indeed that, ip(x) being a smooth function of 
x G M d , and <p(x) being a smooth phase such that V x (p(x) = iff x = 0, and such that 
the Hessian at a; = 0, x </?(0), is invertible, then the stationary phase theorem states the 
asymptotics (see [Ho]) 

r ' (r\ ■ \d/2 

/ dx exp (-(fix) )ib(x) = e% 7777^(0) + 0(e 1+ h . 

V V £ J (det^ xV? (0)) 1/2 

In the case (2.7), we have d = 6. At this stage, we observe that the main contribution in 
(2.8) vanishes, due to the sum over o = ±1. Hence we recover, 

li=0(Ne 3 -^) =0(^/2). 



For Z 2 w e have: 

J 2 = i= V Y * f*dn f ^$(h)e*^- x '^- v '^-^) 

f o, , + ° h ah (2.9) 

/■ (Xi - vtf, . . . ,X r - V r t + — Ti, . . . ,X S - V s t —T 1 ,...,Xj-Vjt, 

J z z 

oh oh . 

V 1, • • • , V r W s + —,..., Vj). 

Now, we compute the limiting behavior of X 2 in a twea/c sense. Hence we test formula (2.9) 
against a smooth function i^(xi, ■ ■ ■ , Xj, v\, ■ ■ • , Vj). We obtain, using the stationary phase 
theorem in the variables x r and h (one could use the variables x s and h as well) 

J dx\ . . .dxjdvi . . . dvj J 2 ip(xi, ... , Xj, vi, . . . , Vj) = 0(^72) = O^e 5 / 2 ^ . 

This again assumes some smoothness and decay assumptions on 
The term X3 can be treated in the same way. We write, 



N-j 

e 

l<r<s<j I 



E E / ^ r dt 1 S{t-T 1 )T? ta S{n-t 1 )C! tj+1 S(t 1 )fi +1 . (2.10) 

r<s<ii=l Jo J ° 



Here, as in formulae (2.7) and (2.9) defining Z\ and X 2 , the T term in (2.10) induces an 
integration over a variable which we call k\ G M 3 (playing the role of h in formula (1.11)), 
and the C term induces an integration over variables which we call hi G M 3 (playing the 
role of h in formula (1.17)), as well as Xj+i G M 3 and Vj+i G 1R 3 . Now, for each value of r, 
s, and the corresponding term in formula (2.10) is tested against a smooth function ip, 
as we did for X 2 . The result involves an integral over all variables x\, . . . , Xj+i, i>i, • • • , 
Vj+i, hi and fci. As in (2.9), we now apply the stationary phase theorem. More precisely, 
if r £, we use the stationary phase in the variables k\ and x r (to handle the T term) , as 
well as hi and a^+i (to handle the C term). In the case r = £, we rather use the variables 
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k\ and x s for the T term, as well as hi and Xj + i for the C term. We refer to the treatment 
of X4 below for details. Using this approach, we eventually obtain, 

'N-j 



x 3 = o(^. 6 )=o( £ 2 ), 



in a weak sense (i.e. in the sense of distributions in xi, . . . , Xj, vi, . . . , Vj). This again 
requires some smoothness and decay properties for the function as well as 0. 

To treat the last term J 4 we have to distinguish various cases: 

1) {r,j + i}n{e,s} = $, 

2) r = £ and s + 1, 

3) r ^£,s = j + l, 

4) r = £, s = j + 1, 

and treat them separately. Item (4) obviously corresponds to a collision/recollision event, 
and this is the dominant term that we want to keep in the limit, while replacing the true 
expansion (1.23) by the reduced series (1.25). 

As regards case 1), we write, 

N-j , f* f 1 ' f f dtn f dkt 



jr,£,s _ a±a[ / dt\ j dr\ / dxj+idvj+i 



e 'Jo l J (2nrJ (2yr)3 

<Tl ,CT 1 =±1 

0(/ii)0(fci)e i ^-(^-*' +1 - Wr ^^ 

a' x h\ a[ki &\h\ (0 11) 

Xi - vgt H — Ti, . . . , x fl - v a * — Ti, . . . ,x r - v r t H — *!, . . . , v • ^ 



- *i! v i» • • • » ^ > • • • » «a H ^— , • • • , "r — , 

As before, we test formula (2.11) against a smooth function ip, and the stationary phase 
theorem in the variables hi, Xj + i as well as ki, xg, readily gives in this case, 



r As = (N_J_ e ^ =0 (e 2 ), 



'iV-j 

e 

in the sense of distributions in xi, . . . , xj, vi, . . . , Vj. 

Case 2) follows along the same lines, except that we now use the stationary phase 
Theorem in the variables h-\ , well as ki, x 8 , in the formula, 

Ms _N-j \- , f* ^ [ ll ,_ f , , f dhi f dki 

-La — 



(7l,(T 1 =±l 

/ i+ i(xi - v x t, xi - v e t H — ti + — ^— ri, . . . ,x s - v s t — n, . . . , 

- Uj+it — h; vi,...,v t — , ...,v s + — — , . . . , v j+1 + 



2 ' 
(2.12) 



10 QUANTUM BOLTZMANN EQUATION 

Case 3) is the same, and we simply need to use stationary phase in the variables h±, 
x r , as well as fci, X£, in a formula analogous to (2.12). All these computation require, as 
usual, some smoothness and decay assumptions on fj +1 and </>. 

Now, case 4) cannot be treated by a direct stationary phase approach, since the space 
variables (the difference xg — a^+i) in the two oscillating factors appearing in this case, 

exp (i— ■ (x£ - x j+1 -v e (t- t 1 )) S ) 

V f J (2.13) 

exp (i— ■ (x e - x j+1 - (ve - v j+1 )(t - ti) - (v e - v j+1 - aihi)(ti - ri))J 

are the same: one cannot decouple the problem into simply applying twice the stationary 
phase theorem (one time for each oscillating factor) as we did before. It turns out that the 
present contribution is, beside X , the only contribution 0(1) which survive in the limit. 

We thus end up this section by carefully analyzing the limit e — > in case (4). The 
analysis performed here allows to understand in a particular case all the arguments needed 
for the general proof given in the next section. We write, 

X 4 " = _ £ a lffl I d tl ^ dn J dx j+1 dv j+1 jj^sj^ys 

<Tl,<J 1 =±l 

0(hi)0(ki)e i ^'( Xe ~ X: > +1 ~ Ve( - t ~ tl ^ e i ^-( x ^- x J+i-(^- l, j+i)(*-*i)-(^-' (; i+i- fJ i /l i)( t i- T i)) 

o-\h\ o\k\ o~\h\ o\k\ 
fj+iix! - vit, ...,x t -v e t+ — g— *i + — T i, • • • , x j+1 - v j+1 t —t x — n; 

v u ...,v t — , v j+1 + + —^-). 

(2.14) 

As explained above, one cannot make use of the oscillations of the product (2.13) by 
simply using the stationary phase in the variables h\, k\, X£ and Xj+\: one has to use 
the oscillations in the velocity variable Vj + i as well. Unfortunately, the factor (ti —T\) in 

(2.13) may vanish, thus killing the oscillation. This is the very reason for the rescaling we 
now perform. We make the following change of variables in (2.14) 

t 1 -r 1 =es 1 , (i.e. Tl =t 1 -es 1 ), 

(2.15) 

£i = (fei + fci)/e- 

The variable s± is the rescaled time between the first "collision" involving particles I and 
j + 1 (occurring at time t\), and the recollision event (occurring at time T\). This gives in 

(2.14) , 

l £ /> j+1 = (N- ,> 3 ^ £ dt i ds ^ j dx j+1 dv j+1 / / 

(7l,CT^=±l 

4>(-h +£^ 1 )0(A;i)e^ 1 '(^" Xj+1 "^ (t " tl) )e- lSlfcl '^-^+ 1 - fJl/ll) /] ) +i(---) • 

(2.16) 
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Now, in order to treat the limit e — > in (2.16), one observes that it is of the form, 

A e = f dxd^dydi] [ ds exp(z£ • x - isrj ■ y) x(x, y, f, rj) , (2.17) 

JR 4d Js=0 

for some smooth function x- We claim that, 

r r+oc 

A £ -> (2tt) 3 / dydri / ds exp(-is V ■ y) X (0, 0, £, rj) , (2.18) 

JR2d Js=0 



with the uniform bound, 



\A e \ < C \\(j r x,yX)( ( ^iPi^'n)\\L 1 (d^dr r ,L^ (dad/3)) 

+ C || (J 7 x,yX)( a ,P,tii r i)\\L 1 (dadf3;L°°(d£dr))) 



(2.19) 



This is merely a variant of the stationary phase Theorem, and it relies on the simple two 
identities, 



dxd£dydr) e ix ^ lS7] - y x(x, V, £, v) = J ^(V^x) (-£, s?7, £, ??) 



(2.20) 



giving, in particular, the absolute convergence of the ds integral in (2.18). This model 
computation proves that the integral (2.16) concentrates asymptotically on the set £i = 0, 
Xj+i = xg + vg(t — ti): particles I and j + 1 eventually collide at the same point. The 
integral over si also becomes an integral over the whole set 1R + as e — > 0, and all factors 
ti = t\ — es\ tend to simply become t\ = t\\ the two collision/recollision events tend to 
take place simultaneously. Note the crucial fact that the concentration in the variables Xj+i 
and £i, stemming from the first oscillating factor, happens independently of the variables 
s\, ki, and Vj+i of the second oscillating factor. Hence we recover the asymptotics, 



a 1 ,a' 1 —±l 



|(^(/ci)| 2 exp ( - isifci • - Uj+i + crifci)) 



J j+1 {xi -vit,...,x e -v e t ti,..., Xj+i -Uj+i*H * 



(2.21) 



2 

(cri - o-i)fei (cri - o-[)ki 
V!,...,V t + , . . . , v j+1 ). 

There remains to use the fact that, 

f oo 



1- 



Re J ds\ exp ^ — is\k\ ■ (vg — Vj+i + o~iki)j = —n5 {k\ ■ (vg — Vj+i + a\ki) j . (2.22) 

Using formula (2.22) together with (2.21), and explicitly performing the sum over o~i, it 
turns out that the term |</>(...)| 2 exp(...) in (2.21) becomes merely 

|^(/ci)| 2 5^/ci • {vg - Vj+i + crifci)) , 
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as we prove in the next section. This factor gives the natural cross-section in the limiting 
Boltzmann equation: it gives the Fermi Golden Rule together with the natural conservation 
of kinetic energy. 

As it is clear, the analysis performed in the present paper heavily relies on the repeated use 
of the stationary phase Theorem, together with asymptotic statements of the form (2.17), 
(2.18), and (2.19). Eventually, everything boils down to checking that various oscillating 
factors involve various independent variables, so as to be able to use that exp(ix 2 /e), 
respectively exp(zsx 2 ), has size e d l 2 , respectively s~ d ^ 2 , in the sense of distributions in 
x G M, d . This kind of observation is standard in the field, and we may quote [Spl], [EY1], 
[EY2], [Cal], [Ca2] for a systematic use. We also redirect the reader to the next section 
for precise statements in the general case. 

3. Convergence. 

The main result of this paper is summarized in the following theorem: 
Theorem 3.1. Suppose 4> G L X (IR 3 ) n L°°(IR 3 ). Assume f is such that 

/ \f (h,k)\dhdk + / dh sup \f (h,k)\ < +oo. 
J J fceR 3 

Then for all j > 1 and for t < to, with to sufficiently small depending on fo and cf>, 

lim fj(t) = f^(t) 

pointwise in Xj,Vj,t. The function fit), t G [0, to) is the solution of the "classical" 
Boltzmann equation 

(d t + v-V x )f = Q(f,f) 
Q(fJ) = J Jdv 1 duB(u J ,\v-v 1 \)(f'f[-ff 1 ) 

with B given by 

B(u, w) = — =\uj ■ w\ \4>(uj (ui ■ w))\ 2 . 

Remarks 

The bounds f° in L 1 (IR d ;L 1 fl L°°(IR d )) in the statement of the Theorem directly stems 
from estimates of the type (2.19) above. 

The assumption cf> G L X (R 3 ) fl L°°(R 3 ) is implied by the stronger </> G Cq(1R 3 ), whenever 
k > 3. 

Finally, the above theorem holds in any space dimension d > 2, as it is clear from the 
proof. 

Proof : For sake of shortness we write explicitly the proof for j = 1. The extension to 
arbitrary j is straightforward, as well as the factorization property of the limit, which 
follows from Lanford's classical argument (see [L] and [CIP]). 
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For j = 1 we can write ff = ff as follows with £\ = 1 

, ^ (N-l)(N-2)---{N-n-l) A ^ /•* , , 

n=0 £ £ 2 =1 C^l 70 " / ° 

/ dt 2 / dr 2 • • • / / dr n S(t - *i)Cf 1>2 S(*i - n)T/ li2 

Jo JO JO Jo 

5(ri - t 2 )Cf 2y3 S(t 2 - r 2 )T| 2)2 . . . S(r n _i - *n)C f !„,„+i5(t„ - T n )T/ njn+1 5(r n )/° +1 . 

(3.1) 

The solution to the Boltzmann equation to be compared with the expansion (3.1) is 



2 n t /.*„_! 

n>0t 2 = l c=i 



(3.2) 



iS(ti — t 2 )CV 2 ,3 . . . S(t n -i — t n )Ci n!n+ iS(t n )f n+1 , 

with 

(C i j +1 fj +1 )(x 1 ,v 1 ,...,x j ,v j ) = dv j+1 du B(u,wt) 

JR3 JS 2 

(3 3) 

. .., a;*, ...,xj, xf,V!, vt-u(u -wt), . .., Uj, Vj+i + u>(w • W£)) 
- /j+i(a;i, a*, ...,Xj, X£-,vi, ...,v £ ,...,Vj, v j+1 ))], 

where § 2 is the unit sphere in IR 3 , = U£ — Vj+i, and -B(cu, is the quantum cross section 
computed in the Born approximation. 

Therefore both f[ and fi are expressed as sums of integrals of f^+i computed in suitable 
points which we construct in the following way: for a fixed n we call graph of order n 
any sequence £i,...,£ n with £j e {1,2, and denote by Yl' e 1 e n ^ e sum on a ^ 

the possible graphs. Therefore upon interverting the t and r integrals in (3.1) for later 
convenience, 



. w (N — l)(N — 2)---(N — n—l). 

tf(x 1 ,v 1 ,t) = J £ E " 7^ -(-*) 



2n 



n>0 £!,...,£„ 
"t ft 1 rt n —i rti pt2 



/ efti / dt 2 - •• dt n dr\ \ dr 2 - • ■ / dr n -i / dr n 

Jo Jo Jo Jt 2 Jt 3 Jt n Jo 

E-A- , f dh\ f dh n f dh f dk n A - , . ln . 
, II J w"J Wf J Wr"'J wn* 1 *' 1 

a 1 ,...,a n ,a 1 ,...,a n =±l J — i J — i 

y dx 2 ---J dx n+1 J dv T -- J dv n+1 f^ +1 (y 1 ,...,y n+1 ;u 1 ,...,u n+1 ) 

JL ft. fe- 
ll exp [i-^ • (tj-) - yj+iitj)) (ye, fa) - y,-+ifa))] ■ 

i=i 

(3-4) 

In (3.4) Vj(s), Uj(s) are the backward trajectory of the particle j '< > 1 and its velocity. The 
particle j is born at time at Xj with velocity and the trajectory is computed for 
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s G [0, tj-i], (as the reader may easily check), by the following formulae: we set to = = t, 
t n+ i = and 

Uj +1 (tj) = v j+1 , 

h h- 

U j + 1 fa )=Uj+l (tj ) + (Tjy U tj (Tj)= U £j (tj ) - Oj -j- , 

k- k- 
u j+1 (t j+1 ) = u j+1 (tj) + a'j-l ui.(t j+1 ) = u e .fa) - a'j-l, 

u r (tj+i) = u r fa) = u r {tj) if r =^£j . 

Also, the values of the Uj : s are extended outside the "collision" times given by the £, 's and 
Tj's, as follows, 

Uj( S )=\ U ^ r) , *T <S<t r , 

\ Uj(t r+1 ) ift r+ i<s<r r . v ' 

The trajectories yj themselves are now simply defined by, 

yj (s) = Xj- / Uj(r)dr. (3.7) 

J s 

With this construction Uj(s) is right-continuous, i.e. at the times t r and r r , the velocities 
are the outgoing ones. We have also set 

yj = yj (0), Uj = Uj(0). (3.8) 

We denote by T e (t l7 . . . , t n ; . . . , i n ) the contribution to the n-th order term of the above 
expansion due to the graph {£i, . . . ,£ n }: 

s2n (N-l)(N-2)---(N-n-l) * 



r e (t 1 ,... ) t B ;f 1) ..,g = (-ir- ^— ^ / d Tl ... 

"tn-1 ft n n 



ptn-1 Pt n " 

/ dr n -! / dr n ^ Ti a J a j 

Jtn J ° a l5 ...,a„,ai,...,<=±l^l 

J dx 2 ---J dx n+1 J dv 2 ---J dv n+1 f° +1 (y 1 ,...,y n+1 ;u 1 ,...,u n+1 ) 

JL h- k- 

[[ exp ■ (y t . (tj) - y j+ i(tj)) + i-f • (vt 3 fa) - y j+1 fa))] , 

3=1 

so that 

i ft rtji-i 

f!(t) = J2J2 / dtl "- dt n T e (t 1 ,...,t n ;£ 1 ,...,e n ). (3.10) 

Analogously, we write 

/i(*) = Z) /^i - / dt n T(t u ...,t n ;e u ...,£ n ), (3.11) 



n>0£ 



1 i ■ ■ ■ j^n 



o Jo 
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with 

T(t 1 ,...,t n ;£ 1 ,...,£ n )= ^2 dv 2 --- dv n+1 du 2 - ■ ■ du n+1 

(3-12) 

( II a o B ^h wj))fi+i(yi, • • • , Vn+i^f, . . . , 

3=1 

where Wj = u\ (tj) — u^ l +1 (tj), and the classical trajectories yf(s) are computed (as the 
reader may easily check), by the following formulae: we set to = t, £ n +i = 0, and 



uf +1 {tj) = v j+1 , 

uf +1 (t j+1 ) = uf +1 (tj) + ^p-{Uj ■ (uf (tj) - Uj+^tj))^ 

(3.13) 

u?(t j+1 ) = U ?(t j ) if r±li. 

Also, the values of the it^'s are extended outside the "collision" times given by the tj% as 
follows, 

uf(s) = uf(t r+1 ) if t r+1 <s<t r . (3.14) 



The trajectories yf themselves are now simply defined by, 

yf(s) = x j - I 3 1 uf(r)dr. (3.15) 

J s 



Note that, because of the factors (1 + oy)/2, when a r = — 1 the velocity is unchanged at 
time t r , so producing a loss term, while when a r = 1, the velocity changes by a term such 
that the conservation of energy is ensured. Finally, 

yf = yf(0), uf = uf(0). (3.16) 

As stated in Theorem 3.1, we want to prove that for all (x±, v±), 

fi N (t)^h(t) (3.17) 

for t sufficiently small, assuming that: 

1) ( ^GL 1 (M 3 )nL 0O (R 3 ) , 

2) The following norms are finite 



Ni(f°) = ii/°iu 1 (r»xr"), mf) = / desup|/°(e, 

J k 



k)\ 



This is obtained from two steps: The first step is to prove that T £ if uniformly bounded 
(Proposition 3.2), the second is to prove Eq. (3.17) via the dominated convergence theorem, 
ensured by Proposition 3.3. 



16 QUANTUM BOLTZMANN EQUATION 

Proposition 3.2: (Uniform bound). There is C > 0, only depending on (p as in 
Theorem 3.1, such that 

\T s {t u ...,*„; ... A)| < C n (^iCf°) + N 2 (f°)r ■ (3.18) 

Proof : The proof relies essentially on (a systematic use of) the change of variables (2.15) 
and the bound (2.19) used to treat the term J^ej+i Q f ^ e previous section. 
First we make the following change of variables 

e hj + ki h^lL. (3.19) 

£ £ 

The expression of T s (ti, . . . , t n ; £±, . . . , £ n ) thus becomes: 
T e (t l7 . . . , t n ; £ u ..., i n ) = {-i) 2n (N - 1)(N - 2) • ■ • (N - n - 1) e 3n 

'(ti-t 2 )/e r(t2-t 3 )/s r(t n -i-t n )/e rt n /e 



/ dsi I ds 2 ■■■ ds n -i I ds n ^ jf[ 
Jo Jo Jo Jo n , 

cri,...,cr n ,a 1 ,...,(T n —±l j — i 

J dx 2 ---J dx n+1 J dv 2 ---J dv n+1 fn +1 (y 1 ,...,y n+1 ;u 1 ,...,u n+1 ) 

n 

JJexp [it, ■ (yi^tj) - yj+iitj))] 



3=1 

n 



[[ exp [ - i^- ■ [( ye . (tj) - y j+ i(tj)) - (y tj (tj - esj) - y j+ i(tj - esj))]] . 

3=1 

(3.20) 

Also, for fixed x\, vi, and fixed £ r , h r , s r , t r , cr r , a' r (r = 1, . . . , n), we define the integral 

/n 
dXdV Yl <j>{kj)(f>(-kj + ££,j)fn+i{yi, Vn+l] Ml, ... , u n+1 ) 
3 = 1 

n 

Y[ exp fa ■ (y e . (tj) - Vj+iitj))] (3.21) 



3=1 

n 



l[exp[-i-r- [{ytj (tj) - yj+i(tj)) ~ (^ fa ~ £S 3) ~ yj+ifa ~ • 

3=1 

Here we use the shorthand notation, 

X = (x 2 , • • • , x„+i) , V = (v2, • • • , v n+1 ) , E = (£i, . . . , f n ) , K = (fei, . . . , k n ) . 

In order to bound T e (or 7), we now follow the same lines as in the treatment of X^'- 7 " 1 " 1 : 
we use appropriate changes of variables so as to come up with factors of the form (2.17). 
Then we treat them as in (2.18), (2.19), (2.20). The description of the relevant changes 
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of variables is the reason for the discussion made in equations (3.22) to (3.33) below. The 
resulting Proposition 3.3 below is analogous to the bound (2.19), and Proposition 3.4 states 
a generalization of the asymptotics (2.18). 

Let us come to the details. We note that the mapping 

(x 2 , i>2, • • • , x n+1 , v n+1 ) — ► (y 2 , u 2 , ■ ■ ■ , y n +i, u„+i) (3.22) 
is one-to-one with unitary Jacobian. Moreover, from Eq.s (3. 7), (3. 8), we have, 

Vj( s ) - (Vj + u j s ) = [ ( u j( T ) ~ u j) dr , ( 3 - 23 ) 

Jo 

so that one readily obtains 

Vis (tj) ~ Vj+iitj) = Vti - Vj+i + ( u £j - u J+i)tj + l), (3-24) 

yii - Vj+iitj) ~ iVtj ( f j - £s j) - Vj+iitj ~ £s j)) = ^j(u £j - u j+1 ) + 7? . (3.25) 

Here 7* , i = 1, 2 do not depend on y r , u r for r = 2, . . . , n + 1. The terms 7^ actually are 
linear in the variables e£ r and h r , for the values r = j,j + 1, . . . , n. This gives in (3.21), 

/n 
dYdU ( J] HkjM-kj + e^^/ftyi, «i)/°(Y, 17) 

(3.26) 

exp [i^j- • (y<>. - y^+i)] exp [i(u^ - u j+1 ) ■ (sjkj + tjQ] , 



where 



r(S,K):=^( 7 j-^- £ - 1 7 J 2 -fc J ), (3.27) 

i=i 

and we use the shorthand notation 

^ = (2/2, • • • , 2/n+l) , *7 = (u 2 , • • • , «n+l) • 

Now, as in (2.20), we explicitly compute the dYdU integral in (3.26). It is crucial at this 
stage that (yi, ui) does not depend upon Y - nor U, and that the phase T/s does not depend 
on Y nor t/ neither. Also, we need the following observation (putting the variables y\ and 
u\ apart): 

n n 

E& • (Wi - w+i) = + (E^^) -fi > ( 3 - 28 ) 

n n 

£fc ' K = ( T5 >^) + (E^^) -«i > (3-29) 

i=i i=i 

n n 
E S 3 k 3 ' ( U 4 - U J+0 = ( S K ' AU ) + ( E S *iJ S 3 k j) ' Ul ' (3.30) 
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where A, T, S are n x n matrices, whose elements are given by: 

A r , s = —S r , s + 8e r ,s+i , (3.31) 

T = diag(ti, . . . ,t n ), S = diag(si, . . . , s n ) , , (3.32) 

and ( • , • ) denotes the inner product in R 3n . The matrix A is upper triangular, with —1 
coefficients on the diagonal. With these notations we write, 

/n _ 
dYdU (j[4>(ki)H-kj + tjj) j^'VfliyuuJf^U) 

3=1 

exp [i(H, AY)] exp [ - i(S K, AU) + i(T S, AU)] . 

Here we have set, 

n 

f(E,K) = r(H, K) + J2 6 hA^3 ■ (Vi " 'jkj ■ «i]. (3-33) 

i=i 

The dYdU integral is now easily computed, 

J = e^/?(ift,ui)(nfe*-*i + ^)) iJ(-A T S,A T ^-A T TS) . (3.34) 

Armed with expression (3.34) for /, we are ready to give uniform bounds on this term, as 
well as convergence results for T e . 

Proposition 3.3. Let us define the function 

n 

g(E, K) := | ( J] <fe)) iJ(-A T S, A T SK - A T T S) , 

i=i 

which is independent of e (the dependence of g upon the Sj 's and tj 's is not made explicit). 
Then, the following estimate holds true, 

/ dzdKg^^Kc^mn+Niinrji-——, ^ 

jR6 n j = l ( 1 + S 3> 

for some constant C > depending on <f>. In particular, I defined by (3.34) satisfies, 

r n 1 

/ 1 d~dK <cHi/°n^(M3 ) (iv 1 (/ o )+iv 2 (/ o )r n 77-7^3 • 

JK. 6n - =1 l 1 "T s j) 

Hence I is uniformly (in e) integrable with respect to the variables s± > 0, . . . , s n > 0. 
Remark 

If the space dimension is d, the factors (1 + Sj)~ 3 above become (1 + Sj)~ d . 
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Proof : 

The following first bound obviously follows from (3.34), together with the fact that 
|detA| = 1: 

J dZdK g(E,K) < C \\4>\\l HR3) \\f%\\ L i mL o° m ) ■ 

On the other hand, to recover decay estimates as the s-variables grow, one makes the 
change of variables K — > S~ X K in (3.34) (valid if the s^'s are, say ^0), and we get, 

/ dZdK g(E, K)<(f[ sf) fdZdK(f[\4>(^) |) \f° n (-A T E, A T K - A T T S) | 

J j=l J j=l 3 

n 

<C (JJ«~ 3 ) ||0||2-(K3) llSllLi(]Rd;Li(]K d )) 1 
J = l 

All this proves the proposition. 



Using the Proposition 3.3 we conclude the proof of the uniform bound (3.18): the factor 

(N - 1)(N - 2) • • • (N - n - l)e 3n 

in (3.20) is bounded by c n , and the number of terms in J2 a cr , is 2 2n . The bound (3.18) 
implies the absolute uniform convergence of the series (3.4) for sufficiently small times. 
Indeed, the number of graphs {£\, . . . ,£ n } is n\, while the time ordered integral gives a 
factor t n /n\. Therefore we have 

oo 

\fi(t;x 1 ,v 1 )\<Y,[C(N 1 (f ) + N 2 (f°))tr , (3.36) 

n=0 

converging for C(AT 1 (/°) + N 2 (f°))t < 1. 

□ 



Proposition 3.4 (Term by term convergence). For all (xi,vi), and any choice of 
t > ti > ■ ■ ■ > t n > and {£x, . . . ,£ n } with £ i G {1, 2, . . 

\imT £ (t 1 ,...,t n -,£ 1 ,...,£ n ) = T(t 1 ,...,t n -,£ 1 ,...,£ n ) . (3.37) 

e— »0 



Proof : We want to pass to the limit e — * in formula (3.20) giving the value of T e . Now, 
equation (3.20) expresses T e as the integral over s\, . . . , s n , together with S and K, of 
the quantity / used above. Hence Proposition 3.3 makes it possible to use the dominated 
convergence theorem, and we can safely interchange the lim e ^o with the integration on 
(si, . . . , s n ), and H, K in (3.20). 

Coming back to formula (3.21) expressing / as an integral over variables X and V, we may 
use that the function /^ +1 (Y, U) (seen as a function of X and V) is a fixed (independent 
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of e) function in L 1 (M 6n ). Hence we can again use the dominated convergence Theorem 
to interchange lim e ^o with the integration on X and V in (3.21). 

From Eq.s (3.5)-(3.7) and using hj = —kj + s£j, we have that yj(s) — > Vj(s), and, for 
s 7^t r , Uj(s) — > Uj(s), where %(s) are defined in a similar way of u C j(s), substituting the 
recursive relation (3.13) with 



u j+l (t j+1 ) = u J+1 (tj) + kj 

u t {t j+1 ) = u t (tj) kj . 



(3.38) 



The trajectories are given by 



3-1 



Vj( s ) = x j- / u A T ) dr ■ ( 3 - 39 ) 



Moreover, from (3.5) and the position hj = —kj + e^j 

ytj(tj) ~ Vi^tj - esj) (_ a 3 



SjU tj (I , - esj) -> Sj (u tj (tj) + Y k J 



y j+1 (tj) - y j+1 (tj - esj) _ ( 
Thus, from (3.20), we obtain 



limT e (£i,...,t n ;£i,...,£ n ) = 

e— »0 



d,S r 



ai,...,a n ,a 1 ,...,a' n —±l J — l 

I dX I 



dS JL . _ (3.40) 

II l^( /c i)| 2 /n+l(2/l,---,Wn+i;'Ul,---,yn+l) 



(2yr) 3n 

J=l 

n 

n ex P * fe • (ye* - - ^ ■ [(% + /c J a j) r j]] ' 

where iDj = U£.(tj) — Uj+\(tj). As indicated at the beginning of the proof, we have now 
come up with factors of the form (2.17). Hence, as in the previous section, there only 
remains to take care of the gE and ds\ . . .ds n integration as we did in (2.18). 

The integration on S produces a delta function concentrating on the set {X \ Xj+i = 
ye (tj),j = l,...,n}, meaning that the particle j + 1 is born at time tj exactly where the 
ancestor is at time tj. Let {yj(t),j = 1, . . . ,n + 1} be the trajectories y~j(t) satisfying this 
additional condition. Since such trajectories do not depend on the s's, we conclude that 

]imT e (t 1 ,...,t n -,£ 1 ,...,e n ) = 

r Air r 



h) 2 " e n'x/^/ # ni«« ..... 

n f-\-oo 

fn+i(yi,---,y n +i;ui, u n+1 ) T [ / cisexp [ - iskj ■ [{wj + kjaj)}] , 

,=1 J ° 
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where the integrals on the s's make sense as distributions (see (2.20) and the bound (2.19)). 

There remains to explicitly compute the integral over s and the sum over the cr's, so as 
to identify T(ii, . . . , t n ; £i, . . . ,£ n ) in (3.41). Using formula (2.22) together with the fact 
that the expression for the limit is invariant under the transformation 

s ► s , kj > kj , <j j ► <jj , <7 j > <j j , 

we can substitute 

p + OO 

I dsexp [ — iskj ■ [(wj + kjaj)]\ with its real part, i.e. n5 (kj ■ (wj + kjcrj)) . 
Jo 

Doing the change of variables rfj = —ajkj, and cr'- = —OjOj, the recursive relations (3.38) 
become ^ 

u j+1 (t j+1 ) = u j+1 (tj) + —^-Vj, 
u e .(t j+1 ) = u tj [tj) - ] ~y ± Vj ■ 
Then the sum on aj is free and Eq. (3.41) becomes 



(3.42) 



Bmr-(t 1 ,...,t B ;/ 1 ,...,£ B )= ^ f[ a d J J dV f[ |0(,,) |' 

/n+l(yi, • • • , 2/n+i; Ml, ... , U n+1 )5 (rfj ■ (Wj ~ Tfj)) . 

We perform partially the integration in 77^ in the following way: let 7 be a smooth function, 
we have, going to polar coordinates, 

//> r+00 
drn(ri)6(ri-(w-Ti))= du d\\ 2 ^(\uj)5(\uj ■ w - A 2 ) 

= / duj (uj ■ w) ^({uj ■ w)uj)\{uj ■ w > 0} 

If 

= — I du \u) ■ w\ 7 ((a; • w)u>) . 

Inserting this identity in (3.42), we recover the classical trajectories defined in Eq.s (3.13)- 
(3.15). Doing the integrations in Xj = \rjj\, we identify the limit in Eq. (3.43) with Eq. 
(3.12), where 

2 



4>((oo • w)ui) 



We remark that the cross-section we have found in the limiting Boltzmann equation is that 
given by the Born approximation for the Quantum Scattering problem. This is known as 
Fermi's Golden Rule (see [AM], [RV], [Bo], [CTDL], [Co] ...). For the low-density limit the 
situation is different. Now all terms in the perturbative expansion of Si nt (t) play the same 
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role: We thus expect that the same result holds as Theorem 3.1, but with a cross section 
given by the full scattering matrix S associated to the potential 0. We mention in passing 
the following well-known fact [RS]: S admits a power series expansion in the potential 0, 
called the Born series expansion, whose first term is precisely given by the Fermi Golden 
Rule, 

s(k) = \m\ 2 +o^) , 

roughly. The weak coupling regime may thus be seen (technically) as some first order 
approximation of the low density regime. We refer to [Nil], [Ni2] for a kinetic interpretation 
of the scattering matrix in a semi-classical regime. We also refer to [Cal], [Ca2], for the 
derivation of a linear Boltzmann equation involving the Born series expansion in a low 
density regime. We finally quote [Co] for a physical discussion of related questions. 

Some comments are in order. 

We underline once more that while in the low-density regime classical and quantum 
systems evolve similarly according to the Boltzmann equation, the situation changes dras- 
tically in the weak-coupling limit. Here contrary to the behavior outlined in the present 
paper for quantum systems, classical systems of particles are expected to satisfy the Lan- 
dau equation which a diffusive character (see [Sp2]). Unfortunately no rigorous result is 
known in this direction. The situation is better understood for linear problems. Obviously 
the same scalings can be also considered for a test particle moving in a random distribu- 
tion of obstacles. In this context it is possible to derive a linear transport equation for a 
classical particle in the low-density regime (see [G], [BBS], the review paper [Sp3], or the 
textbook [Sp2]) or the linear Landau equation in the weak coupling limit (see [DGL]). A 
linear Boltzmann equation can be derived in the weak-coupling limit either for short times 
(see Refs [Spl],[La], [HLW] ) or globally in time (see [EY1], [EY2] ) . The one-dimensional 
case is somehow pathological (see [EPT]). For the case of non random scatterers negative 
results are available [CP1], [CP2] (see also [BGW] in the classical context). 

Similar considerations in the case of an atom coupled to a gas have been developed in 
[Du]. In addition, notice that the problem of the wave motion in a random medium is 
of interest for the applications as shown in Ref.s [KPR] We finally mention [PV] for the 
analysis of a weak coupling regime when the obstacles are temporally random (and the 
underlying process is at once (almost) Markovian) . 

At the physical level, the question of passing from the Schrodinger equations to (linear 
or non-linear) Boltzmann equations is an old problem. We may quote [Pa], [KL1], [KL2], 
[Ku], and [VH1], [VH2], [Zw], as well as the textbooks quoted before. 

When this paper was finished we were informed by H.T. Yau that arguments similar to 
those of the present paper have been independently developed by him and his group in a 
paper in preparation. 
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